Estimation of rock Fe content based on hyperspectral 
indices 
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Abstract: Information on the Fe content of bare rocks is needed for implementing geochemical processes 
and identifying mines. However, the influence of Fe content on the spectra of bare rocks has not been 
thoroughly analyzed in previous studies. The Saur Mountain region within the Hoboksar of the Russell 
Hill depression was selected as the study area. Specifically, we analyzed six hyperspectral indices related to 
rock Fe content based on laboratory measurements (Dataset I) and field measurements (Dataset II). In situ 
field measurements were acquired to verify the laboratory measurements. Fe content of the rock samples 
collected from different fresh and weathered rock surfaces were divided into six levels to reveal the spatial 


distributions of Fe content of these samples. In addition, we clearly displayed wavelengths with obvious 
characteristics by analyzing the spectra of these samples. The results of this work indicated that Fe 
content estimation models based on the fresh rock surface measurements in the laboratory can be applied 
to in situ field or satellite-based measurements of Fe content of the weathered rock surfaces. It is not the 
best way to use only the single wavelengths reflectance at all absorption wavelengths or the depth of these 
absorption features to estimate Fe content. Based on sample data analysis, the comparison with other 
indices revealed that the performance of the modified normalized difference index is the best indicator 
for estimating rock Fe content, with R? values of 0.45 and 0.40 corresponding to datasets I and II, 
respectively. Hence, the modified normalized difference index (the wavelengths of 2220, 2290, and 2370 
nm) identified in this study could contribute considerably to improve the identification accuracy of rock 
Fe content in the bare rock areas. The method proposed in this study can obviously provide an efficient 
solution for large-scale rock Fe content measurements in the field. 
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1 Introduction 


Iron (Fe) minerals in rocks play a pivotal role in geochemical investigations; for this reason, it 
is necessary to monitor Fe content of rocks in extensive areas (Baugh et al., 1998; Velasco et 
al., 
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2005; Li et al., 2009; Liu et al., 2010; Simandl and Paradis, 2018). Conventionally, local 
measurements are an important means for determining Fe content of rocks, and extrapolating the 
measured value from each sampling point is helpful for deriving maps of rock Fe content (Liu et 
al., 2016). Application of remote sensing technology to obtain the feature information of rocks 
has the advantages of wide coverage and fast speed (Hao et al., 2019; Zeraatpisheh et al., 2019). 
Furthermore, some scholars have emphasized that the optical domain can be used for such 
researches (Galvao et al., 2008; Middleton et al., 2011; Nair and Mathew, 2012; Metelka et al., 
2015; Qian et al., 2019). Accordingly, the retrieval of hyperspectral images provides a new 
approach for extracting surface information (Clark and Roush, 1984; van der Meer, 2018). 
However, classic remote sensing instruments (such as the ASTER and Landsat platforms) cannot 
effectively retrieve surface characteristics of rocks because it is difficult to analyze their 
diagnostic spectral characteristics within the wavelengths of 100—200 nm bands (Clark and Roush, 
1984; Wu et al., 1997; Liu and Vekerdy, 2003). In contrast, hyperspectral indices, e.g., leaf area 
index (LAI), chlorophyll concentration, and other plant pigments (Rathod et al., 2016), plant 
water content (Pour et al., 2019), soil moisture (Sarathjith et al., 2016), and soil salinity (Jia et al., 
2021), exhibit very good performance at acquiring information on vegetation and soil 
characteristics. Nevertheless, although some scholars have employed hyperspectral technology to 
quantitatively determine the contents of rock elements, few have leveraged this potential well 
(Sracek et al., 2004; Gasmi et al., 2018). 

Rock reflectance has been shown to be affected by Fe content in the wavelengths of 2000—2500 
nm. Although some studies have indicated that the general trend of rock spectral reflectance 
changes with variation in rock Fe content, some other results suggested that rock Fe content can 
be precisely estimated by hyperspectral data (Clénet et al., 2011). However, this theory is not yet 
fully understood with regard to the estimation of bare rock element contents, and this method has 
a strong innovation in remote sensing application. 

Due to the differences of Fe content in different rock types, spectral absorption should be fully 
considered in the estimation of rock Fe content (Clark et al., 2003; Wang et al., 2019). Rocks are 
composed of many different materials containing various kinds of minerals. Hence, the 
reflectivity of rocks depends on the inherent scattering and absorption characteristics of the 
components and arrangement of rocks; consequently, the influence of rock Fe content on the 
reflectance spectrum cannot be accurately described (Chabrillat et al., 2002). This study is one of 
the first to evaluate Fe content of rocks using hyperspectral data. 

The main purposes of this study were: (1) to explore both the variations in rock Fe content in 
the field and the changes in the spectral reflectance at six wavelength levels; (2) to identify 
hyperspectral indices of rock Fe content that are insensitive to other rock components; and (3) to 
develop a method for estimating Fe content in the field using an optimized Fe content 
hyperspectral index. To achieve these goals, we conducted field surveys as experiments, obtained 
spectral data by collecting bare rock samples and analyzing Fe content in the laboratory, and 
measured and analyzed the spectra characteristics of the samples. The results are of great practical 
significance for the determination of rock Fe content. 


2 Materials and methods 


2.1 Study area 


The study area is located in the Saur Mountain region within the Hoboksar of the Russell Hill 
depression (46°25'—46°52'N, 85°00’—86°15'E). This area has a continental climate, with wide 
temperature range and sparse vegetation coverage. There are a large number of exposed rocks at 
the selected sampling points. The study area contains evidence of coal, copper, siderite, magnetite, 
limonite, titanium magnetite, chromite, pyrite, fluorite, clays, and other minerals in a total of 40 
useful deposits, and five diffusion halos are discovered comprising numerous beneficial minerals 
(elements), including zircon, hafnium zircon, chromium, nickel, and cobalt. 
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2.2 Reflectance measurements 


Both weathered surfaces and fresh surfaces were used for reflectance measurements. Reflectance 
spectra were measured in situ for each sample using an Analytical Spectral Device (ASD) FieldSpec 
FR spectrometer (American Analytical Spectral Instruments Company, Tucson, USA), covering 
wavelengths from 350 to 2500 nm at an interval of 1 nm. The key sampling area was determined in 
combination with remote sensing images and geological data. A total of 313 rock (81 samples) 
and mineral (232 samples) samples were collected. The rock reflectance spectra were measured in 
situ within the period from 10:00 to 14:00 (LST) during 3 August—21 September in 2019 under 
cloud-free conditions with a 25° fiber optic. In each measurement, the radiance and reflectance were 
converted using the "reflectance spectrum" function on the panel. These in situ field reflectance data 
were denoted as Dataset II. All spectral measurements were conducted in a completely dark room to 
avoid contamination by other light sources. This laboratory reflectance data were denoted as Dataset 
I. We further divided the laboratory and field sample data according to the degree of weathering. 
Specifically, Reflectance-Laboratory-Old refers to the weathered surface in the laboratory, 
Reflectance-Laboratory-New refers to the fresh surface in the laboratory, Reflectance-Field-Old 
refers to the weathered surface in the field, and Reflectance-Field-New refers to the field fresh 
surface. 

After acquiring the field spectral measurements, the rock samples were transported to the 
laboratory; then, the surface reflectance of each rock sample was again measured using the ASD 
FieldSpec FR spectrometer. A tungsten quartz halogen filament lamp was used to illuminate the 
rock samples from 50 cm away, and the beam was oriented 30° from the vertical. The reflected light 
was collected perpendicular to the sample with a 5° field of view at a distance of 10 cm. We carried 
out measurement experiments in a completely dark laboratory to avoid the effects of other light 
sources. Based on the results of 81 rock samples, we classified Fe content into six levels: 
0.00%-—1.00% (Level 1), 1.00%-—2.00% (Level 2), 2.00%-—3.00% (Level 3), 3.00%—4.00% (Level 
4), 4.00%-5.00% (Level 5), and >5.00% (Level 6). 


2.3 Continuum removal analysis 


To analyze the effects of Fe content on rock reflectance, we clearly displayed the absorption 
characteristics in the rock reflection spectrum by optimizing the continuous measurement spectrum 
(Clark and Roush, 1984). This method was the first to recommend employing a continuum removal 
analysis to distinguish individual absorption features of interest. The continuum is a convex "hull" 
fitted to the straight-line segments at the top of the spectrum. These spectral segments are connected 
to the local spectral maximum and represent the "background" absorption. The other absorption 
features are then superimposed on the background. The continuum can be removed by dividing the 
reflectance value of each point in the absorption feature using the reflectance level along the convex 
hull at the corresponding wavelength. 


2.4 Calculation and determination of the best spectral index 


This research employed six spectral indices to quantify Fe content, as described in Table 1. First, 
according to datasets from the field experiment, we performed correlation analysis to screen the 
combined wavelengths to determine the best wavelength for each index, and then applied Dataset 
II to verify the results of the previous step to determine the optimal index category and the most 
suitable wavelength. The screening step was 5 nm, and the indices were calculated for each 
combination. 

Correspondingly, linear and multiple regression analysis methods can be used to determine the 
relationship between rock Fe content and reflectivity, and thus predict rock Fe content. In this 
study, the model used to estimate Fe content was verified from the aspects of stability and 
prediction capacity. Stability was tested with the adjusted coefficient of determination (R°) while 
prediction capacity was judged with the root-mean-square error (RMSE). R? and RMSE indicate 
the effect of the model. The higher the R? value is, the higher the stability of the model; the 
smaller the RMSE value is, the higher the accuracy of the model. In addition, the relative RMSE 
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(average of the RMSE of Fe content; %) was used to evaluate the accuracy of the regression 


model. 
Table 1 Descriptions of the generic types of indices 
Index Definition Formula Reference 
R Reflectance Pu 
D Reflectance difference index Pr—Pr2 Le Maire et al. (2004) 
SR Simple ratio Pa2/ Pai Le Maire et al. (2008) 
ND Normalized difference index (Pri—Pa2)/(Prrtpr2) Wang et al. (2011) 
mSR Modified simple ratio (Pri-Pa2)/(Prrtpr2-2 paz) Sims and Gamon (2002) 
mND Modified normalized difference EET Sims and Gamon (2002) 


index 


Note: p refers to the reflectance value in the experiment; A1, 42, and (3 represent wavelengths. 
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RMSE = |X (y! -y,} /n, (2) 


j=l 
Relative RMSE (%)=RMSE/y x100% , (3) 


where y' is the estimated value; yis the mean value of y; y; is the independent reference 


measurement; and n is the number of samples. 


3 Results 


3.1 Variation in Fe content of rock samples 


Fe content of 81 rock samples varied greatly, ranging from 0.24% to 8.37%, with an average of 
2.26% and a standard deviation of 1.62%. As mentioned earlier, we divided Fe content into six 
levels: 0.00%-1.00% (Level 1), 1.00%-2.00% (Level 2), 2.00%-3.00% (Level 3), 3.00%—4.00% 
(Level 4), 4.00%-5.00% (Level 5), and >5.00% (Level 6). As can be seen from Figure la, Fe 
content of 81 rock samples can be categorized into three plots. Rock Fe content exhibited broad 
spatial variability, especially in Plot 1. As shown in Figure la, the samples in Plot 1 were 
concentrated, but Fe content of these samples varied from Level 1 to Level 6. In addition, for Plot 
2 and Plot 3, Fe content of the samples was also relatively concentrated. 

Level 2 had the largest number of samples, followed by levels 1 and 3, while levels 4—6 
contained the fewest samples (Fig. 1b). As shown in Figure 1b, all of the rock samples in Plot 3 
belonged to the levels 1-3, except for one sample belonging to Level 5, while a total of 15 
samples in Plot 1 and Plot 2 belonged to the levels 4—6. 


3.2 Spectral properties of rock samples 


For all four reflectance groups (the reflectance of fresh and weathered rock surfaces was 
measured in the laboratory and in the field), the measured average spectra of all 81 rock samples 
displayed similar characteristics, with obvious shoulder peaks at 800 nm wavelength and obvious 
absorption characteristics at wavelengths of 1400, 1900, 2250, and 2350 nm (Fig. 2). The depth 
and width characteristics of the absorption wavelengths at 1900 and 2350 nm were more obvious 
than those of the other two wavelengths (1400 and 2250 nm). The absorption characteristics at 
wavelengths of 2000 and 2050 nm measured in the field were clearly weaker than those measured 
in the laboratory. Moreover, the reflectance values measured in the field were larger than those 
obtained in the laboratory for both fresh and weathered rock surfaces. By comparing the spectra 
between the fresh and weathered rock surfaces, it was found that they exhibited almost identical 
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Fig. 1 Fe content of 81 rock samples. (a), spatial variation of Fe content of 81 rock samples for the three plots; 
(b), number of samples in each level. Plot 1, Plot 2, and Plot 3 represent the spatial aggregation distribution of 81 
rock samples. Level 1, 0.00%-—1.00%; Level 2, 1.00%—2.00%; Level 3, 2.00%-3.00%; Level 4, 3.00%—4.00%, 
Level 5, 4.00%—5.00%; Level 6, >5.00%. 


reflectance data at wavelengths of 350-1000 nm, and the fresh rock surfaces had larger 
reflectance values than the weathered rock surfaces at wavelengths of 1000-2500 nm. 

When the rock Fe content was divided into six levels, the spectra also varied among these 
levels, especially the absorption features at wavelengths of approximately 2250 and 2350 nm 
(Fig. 1). The depths of the absorption features at these two wavelengths increased when Fe 
content increased, especially in the Level 6. After utilizing the continuum removal method, the 
characteristics became more obvious. The amplitude of the absorption feature at the wavelength 
of approximately 2250 nm was small, while that of the absorption feature at the wavelength of 
approximately 2350 nm was much larger. The relative depth difference of the absorption 
characteristics between Level 1 and Level 2 was small near 2350 nm, while the depths of the 
absorption characteristics in the levels 3-5 were more obvious. Further, Level 6 had the greatest 
depth of the absorption characteristics (approximately twice as deep as the other depths) at this 
absorption wavelength. 


3.3 Spectral properties of rock samples 


Table 2 shows the validation of Dataset II based on Dataset I through the reflection wavelength 
with the highest R? and the lowest RMSE. Values of R? and RMSE for the two datasets were 
verified, and the model for predicting Fe content was calibrated using Dataset I (Table 2). 
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Fig. 2 Reflectance properties of rock samples. (a), mean reflectance of the four groups, including Reflectance- 
Laboratory-Old, Reflectance-Laboratory-New, Reflectance-Field-Old, and Reflectance-Field-New; (b), mean 
reflectance of the six Fe content levels; (c), continuum removal reflectance of the six Fe content levels within the 
wavelengths of 2100-2500 nm. 


For the fresh rock surfaces, the reflectance index at 2210 nm exhibited the lowest R (0.05) and 
the highest RMSE (1.64) based on both datasets (Table 2) by comparing the six indices. For 
Dataset I, the performance of the reflectance difference index was moderate (the best wavelengths 
of 2220-2240 nm, R?=0.29, and RMSE=1.47), but for Dataset II, its performance was slightly 
lower (R?=0.25 and RMSE=1.55). Through experimental calculations and analysis, the 
performance of the best normalized difference index at 2250 and 2350 nm was relatively excellent 
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(R°=0.41 and R7=0.38 for datasets I and II, respectively). The performance of the best simple ratio 
index was similar to that of the best normalized difference index, with the same wavelengths and 
only slightly lower accuracy (R?=0.40 and R?=0.37 for datasets I and II, respectively). Modified 
normalized difference index (the wavelengths of 2220, 2290, and 2370 nm) performed best 
among the six indices, with R? values of 0.45 and 0.39 and RMSE values of 1.24 and 1.34 for 
datasets I and I, respectively. The best modified simple ratio index yielded a similar result to the 
best modified normalized difference index, with three wavelengths at 2210, 2250, and 2360 nm, 
although the accuracy was slightly lower (R’=0.42 and RMSE=1.26 for Dataset I and R’=0.38 and 
RMSE=1.37 for Dataset II). 

Measurements of the weathered rock surfaces were consistent with those of the fresh rock 
surfaces for all six indices (Table 2). The best indicator for the weathered rock surfaces was 
modified normalized difference index, with R? value of 0.40 and RMSE value of 1.26 for dataset 
I; however, for Dataset II, its performance was slightly reduced, with R? value of 0.38 and RMSE 
value of 1.37. This index had the same wavelengths for the weathered rock surfaces as that for the 
fresh rock surface measurements, and the accuracy was similar, albeit just slightly lower than the 
wavelengths of the fresh rock surface counterpart. However, it is worth noting that the regression 
coefficient was unstable and its range varied between the fresh and weathered rock surfaces. 

In the experiment of estimating rock Fe content, the determined index was used in the 
exponential regression equation. We used the estimated Fe content and the measured data to 
construct the scatter plots, as shown in Figure 3. It can be seen that most of the data points were 
distributed along the 1:1 line, which suggested that this index had good performance for both field 
and laboratory measurements. 


Table 2 Results of the general types of indices calibrated with the Dataset I and validated to the Dataset II for 
the estimation of rock Fe content 


Beck B Wavelength (nm) Feet bxindex Dataset I Dataset II 
type al V = 23 a b R RMSE e ) R RMSE eae ) 

R 2210 - - 0.43 0.01 0.03 1.64 72.57 0.05 1.63 72.13 
D 2220 2240 - 0.48 0.24 0.30" 1.42 62.84 0.25* 1.54 68.15 
n SR 2220 2250 - 3.64 4.06 0.35" 1.30 57.53 0.32" 43 63.28 
surface ND 2220 2250 - 0.35 6.86 0.36" 1.31 57.97 0.33" 43 63.28 
mSR 2210 2250 2360 0.47 2.20 0.38" 1.27 56.20 0.34* 39 61.51 
mND 2220 2290 2370 0.21 0.79 0.40° 1.26 55.76 0.38" 37 60.62 
R 2210 - - 0.51 0.01 0.05 1.64 72.57 0.06 1.63 72.13 
D 2220 2240 - 0.68 0.27 0.29" 1.47 65.05 0.25° 55 68.59 
Fresh SR 2220 2250 - 4.55 5.07 0.40" 1.30 57.53 0.37" 1.36 60.18 
ee ND 2220 2250 - 0.64 7.24 0.41* 1.30 57.53 0.38* 1,32 58.41 
mSR 2210 2250 2360 0.54 2.16 0.42* 1.26 55.76 0.38" 1.37 60.62 
mND 2220 2290 2370 0.34 0.62 0.45" 1.24 54.87 0.39" 1.34 59.30 

Note: “-”, no data; *, significance level at P<0.05. Dataset I, laboratory measurements; Dataset II, in situ field measurements. R?, 


coefficient of determination; RMSE, root-mean-squared error, Relative RMSE, relative root mean-squared error. 


4 Discussion 


4.1 Spectral properties of rock samples 


Although there are many high-quality optical remote sensing instruments in orbit, few scientific 
investigations were conducted using different wavebands (including visible light, far-near 
infrared, and/or hyperspectral instruments) to accurately estimate the mineral (element) contents 
in rocks, including Fe content. This is because electromagnetic waves cannot completely 
penetrate clouds and vegetation canopy and are further affected by atmospheric factors. The 
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Fig. 3 Comparison between the measured and estimated Fe contents based on the modified normalized 
difference index. (a), fresh rock surfaces with laboratory reflectance; (b), fresh rock surfaces with field 
reflectance; (c), weathered rock surfaces with laboratory reflectance; (d), weathered rock surfaces with field 
reflectance. 


factors influencing the measurements of rock reflectivity include not only Fe content but also 
contents of other elements in rocks, the physical structure of rocks, and the observation conditions 
and methods. Due to these reasons, the observations of Fe content in mixed rock samples are 
biased. In this research, the conditions of the natural arid environment of western China, 
including the measurable bare rocks, sparse vegetation, and relatively cloudless skies, are all 
conducive to calculate Fe content of rocks in the study area. 

The results of rock measurements by previous researchers in different environments have 
shown that Fe ions have major absorption features at wavelengths of 2260 and 2350 nm (van der 
Meer et al., 2014). It has also been found that the depth of these absorption features increases with 
increasing rock Fe content (van der Meer et al., 2014), which is consistent with the results 
reported in this study. Specifically, as shown in Figure 2, when Fe content level increased, the 
depth of the absorption features at wavelengths of 2260 and 2350 nm increased monotonically. 
However, when using only the single-band reflectance at these absorption wavelengths or the 
depth of these absorption features to estimate Fe content, the result was not ideal. This may be 
because the effect of Fe content on rock reflectance is weak, nonlinear, and complicated. 
Furthermore, this effect is influenced by interference due to other elements, such as magnesium 
(Rowan et al., 2004). Therefore, rock Fe content must be identified from hyperspectral indices 
that are insensitive to other factors. 


4.2 Hyperspectral indices to estimate rock Fe content 


Due to differences in applications, some multispectral indices are designed to be insensitive to 
changes in some spectra. For example, the reflectance difference index is not sensitive to the 
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cumulative change in reflectance, while the simple ratio and normalized difference index are not 
sensitive to the proportional effect. Additionally, the modified simple ratio and modified 
normalized difference index are insensitive to both additive and proportional effects (Le Maire et 
al., 2008). Referring to Dataset I, the result indicates that the performance of the modified 
normalized difference index and modified simple ratio is better than that of the simple ratio, 
normalized difference, and reflectance difference index, and the index efficiency considering both 
the additive effect and the proportional effect will be better than the index efficiency considering 
only one effect. The same is true for the Dataset II. The results obviously demonstrate that more 
complex indices, such as the modified simple ratio and modified normalized difference index, are 
better than the other simple indices, which may be because the indices using three wavelengths 
containing more information and thus are more sensitive to other rock properties than those using 
two wavelengths. Therefore, these higher-complexity indices are more appropriate for devising 
more efficient indices. 

Hyper spectral indices are widely used in retrieving vegetation and soil properties (Zeng et al., 
2021). However, few studies have used this approach to estimate the bare rock Fe content. This 
study aims to develop an optimal hyper spectral index for estimating rock Fe content that is 
feasible for both laboratory and in situ field measurements. The best indicator that was identified 
in this research work is the modified normalized difference index, which demonstrates the best 
performance with both laboratory-based measurements and field-based measurements. 
Furthermore, as shown in Table 2 and Figure 3, modified normalized difference index (2220, 
2290, and 2370 nm) performed consistently among the different datasets and rock surfaces (fresh 
and weathered rock surfaces), although there was a small decrease in R? for the Dataset II from 
the weathered rock surfaces. 


4.3 Reflectance of the weathered and fresh rock surfaces measured in the laboratory and 
field 


A large calibration dataset is very important to obtain general indices. Because of cost constraints, 
a restricted research area can have a substantial impact on the experimental design and field 
measurements; as a result, it is impossible to provide a more accurate comprehensive field 
dataset, and only the index method suitable for the current research area can be used. By contrast, 
fitting experimental datasets to actual field-measured spectra can provide an alternative solution 
to find a suitable index. However, to use laboratory datasets instead of field measurements, the 
laboratory-measured rock reflectance spectra must completely match the field-measured spectra. 
Under general conditions, based on the measured spectra, the matching results between the 
laboratory and field spectra are not very satisfactory. The experimental results of this research 
prove the feasibility of the proposed method, which accurately describes the correlation of the 
rock reflectance spectra measured in the laboratory and in the field. Consequently, some 
characteristic wavelengths can be identified with the proposed method (for example, the 
2000-2100 nm wavelengths, and the field-measured reflectance is always larger than the 
laboratory-measured reflectance at all optical wavelengths. Therefore, all regression models based 
on the laboratory datasets may need to be recalibrated before being applied in the field. Thus, 
extending actual spectral data obtained in the field by using spectral data from the laboratory is an 
excellent solution. However, due to the shortcomings of this approach, it is difficult to find 
suitable accurate indicators, so the comprehensive utilization of laboratory and field data can 
verify and calibrate the final results. 

In recent years, the methods employed to accurately and quantitatively determine rock Fe 
content are unable to combine laboratory measurements with in situ field measurements 
(Moumane et al., 2021). In this research, two datasets comprising laboratory measurements and 
field measurements were comprehensively utilized, and the correlations between bare rock Fe 
content and spectral indices were accurately and extensively described and tested. 

Many laboratory measurements use the fresh rock surfaces of rock samples, but in situ field 
measurements and satellite-borne sensors have access only to the weathered rock surfaces when 
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obtaining reflectance information. In this study, we measured the reflectance of experimental 
samples from both weathered and fresh rock surfaces to detect the difference between the two 
types of surfaces. The results indicate that the patterns of the reflectance curves of the two types 
of surfaces were almost the same not only for the laboratory datasets but also for the field 
datasets; the only difference was that the reflectance of the weathered rock surfaces was slightly 
higher than that of the fresh rock surfaces. The best Fe content estimation indices for these two 
surfaces were also the same (Table 2), and the accuracy was just slightly lower for the weathered 
rock surfaces than for the fresh rock surfaces (the best R? values of 0.45 and 0.40 for the fresh and 
weathered rock surfaces, respectively). This small error is to be expected and acceptable, as the 
weathered rock surfaces always have adventitious materials, which will affect the estimation 
accuracy of Fe content. The outcomes of this work indicate that Fe content estimation models 
based on the fresh rock surface measurements in the laboratory can be applied to the field or 
satellite measurements of the weathered rock surfaces. 


5 Conclusions 


This study assessed the use of hyper spectral indices for estimating the bare rock Fe content. The 
longwave-infrared domain (1100-2526 nm) is suitable for predicting rock Fe content, especially 
in two or more characteristic wavelengths, which is somewhat obvious. The best indicator for 
estimating rock Fe content in the field is the modified normalized difference index (2220, 2290, 
and 2370 nm), a commonly used spectral parameter that can accurately and quickly measure Fe 
content in bare rocks. The method proposed in this study can obviously provide an efficient 
solution for large-scale rock Fe content measurements in the field. In order to improve the 
accuracy of measuring rock Fe content, it is necessary to accurately describe the correlation 
between Fe content and reflectivity in rocks and to ascertain how different characteristics of rocks 
and elemental contents in rocks affect their Fe content. 
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